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Crushing candies on the line 
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1 The candy crush game 

In this paper we investigate stability properties of a probabilistic cellular 
automaton. The model is based on the candy crush game (see e.g. nua). 
The idea of the game is that each point (these are the candies) in a rect¬ 
angular grid has a certain color. All candies making a horizontal or vertical 
monochromatic chain disappear simultaneously and candies fall from the top 
to fill in the gaps. If the resulting configuration again contains such chains, 
a sequence of reactions occurs and the player is awarded bonus points. This 
game inspired us to study a somewhat simplified model, where one of the 
main questions is under what conditions infinite sequences of reactions can 
occur. 

2 Model based on the game 

Our model is defined on Z rf with the Euclidean metric. The points in Z rf 
are called sites. We choose n different colors and define the color set C = 
{ci,...,c n }. Furthermore, we choose a stability constant k £ N, k > 2. A 
coloring or configuration is a map rj : Z rf — > C. Given the coloring, we will 
define the notion of stability. Sites aq,... ,x m £ TL d satisfying 

Xi - xt-i = x 2 - xi, ||aq - Xj_i|| = 1 and rj(xi) = rj(x i) 

for i = 2,..., m are said to make a monochromatic chain of length m. Sites 
that are in a monochromatic chain of length greater than or equal to k are 
called unstable. All other sites are stable. Given a configuration r/ we define 
a stability function a v : 7L d —> {0,1} by 


O'lfix ) is stable in r;}- 

A coloring rj is stable if and only if t7 r? (x) = 1 for all x £ Z d . A stable 
configuration in which the color of each site only depends on the parity of 
the sum of its coordinates will be referred to as a chessboard coloring. The 
set of all configurations will be denoted by A. 
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In our model, the dynamics are simpler than in the game. We will use a time 
variable t, taking values in N. The configuration at time t will be denoted 
by rjt , so t]q is the initial configuration. All sites that are unstable at time t 
(that is, sites x for which cr Vt (x) = 0) will be recolored simultaneously and 
independently according to some probability distribution p on C to construct 
qt +1 ■ For i] £ A, we will denote the random configuration that results from 
recoloring the unstable sites by R(rj), so rjt +1 = R{jj)- 

3 The one-dimensional case with two colors 

Let d = 1 and k = 3. Furthermore, let the color set be coded by C = {0,1} 
and choose equal recoloring probabilities, so p = (^, ^). Then we have the 
following result: 

Theorem 1 Choose an unstable initial configuration i] o- Suppose there ex¬ 
ists M £ N such that cr Vo (x) = 1 for all |x| > M. Then 

P p (3to : r /1 is stable for all t > to) = 1. 

Proof. First define the number of unstable sites at time t: 

It = | {x : a Vt (x ) = 0} |, t > 0. 

Note that R is bounded by 2M +4t+l for all t. We will show that linp-j.oo R = 

0 a.s. We define an upper bounds for the probability that an unstable site 
is instable again after k time steps as follows: 

p{ = sup J?eA p (<TRk(„)(x) = 0 I a v (x) = 0), 

p[ n = sup )?eA p ((TRk^ix) = 0 I a v (x - 1) = a v (x) = a v (x + 1) = 0). 

By translation invariance, these upper bound do not depend on x. Note that 
p{ H < Pk- a similar way we define upper bounds for the probability that 
a stable site is instable after k recolorings of the configuration. Since this 
probability highly depends on the distance to instable regions, we condition 
here on stability of a neighborhood of x: 

m m+1 

pf(n,m) = supP(a R k^(x) = 0 | Q <J v (x+i) = 1, Y'' a^x+i) = n+m+1), 

T7GA n 

1 i=—n i——n— 1 

where n, m are non-negative and allowed to take the value oo. First we 
remark that these probabilities are symmetric: pf(n,m) = p^.(m,n). Since 
a stable site can not get instable if a large enough neighborhood is stable, 
the probabilities p)?(n, m) satisfy the following property: 

pf(n,m) = 0 if n,m>2k. (1) 
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Furthermore, if there are at least 2k stable sites at one side, pf (n, m) does 
not depend on the exact number of them: 

pf(ni,m) = pf(n 2 ,m) = pf(oo,m) if ni,n 2 >2k, . , 

pf(n,mi) = pf(n, m 2 ) =pf(n,oo) if mi,m 2 >2k. 

Let 77 € A. A set {x, x + 1,..., x + 5 } is called a bounded stable region of 
size 5 if a v (x) = ... = a r] (x + g) = 1 and cr^(x — 1 ) = cr J7 (x + g + 1 ) = 0 . 
Suppose G is a bounded stable region of size g, then the expected number 
of sites in G that get instable in k steps is bounded from above by 

9 

Y p k^ ~ l ,9 -*)■ ( 3 ) 

i— 1 


This sum is the same for all g > 4fc, since in that case 


9 

2=1 


2k g 

- l,oo) + X Pk(°°i9-i) 

2=1 2=0—2fc+l 

2k 4fc 

i » 4fc -*)+ XI 

2=1 2 = 2 / c +1 

4/c 

Y p k(' 1 ~ - *)> 

2=1 


by the properties © and ([2|). If g is such that a v (x) = 1 and er, 7 (y) = 
0 for all y < x, then the set {y £ Z : y < x} is called a left-unbounded 
stable region. A right-unbounded stable region is defined similarly. For an 
unbounded stable region, the expected number of sites that gets instable is 
at most 


00 2k 4 k 

Yd(i-i,°o) = Y^(i-i,°o)= X pU 4 k ~h°°)- 

i= 1 i =1 i=2k+l 


where we used (pQ) for the second equality. Exploiting (j2]h we arrive at 


00 

YPkd - !'°o) 

2=1 


^ 2k ^4 k 

2^ZPk(. i ~ 1 Ak - i) + - X pf(4fc-*,*- 1 ), 

2=1 2=2fc+l 

^ 4/c 

2 X^(* — i,— *). 

2=1 


Since an instable region consists of at least 3 sites, at least a third of the 
unstable sites has two unstable neighbors. For the same reason, the number 
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of bounded stable regions in rjt is at most It/ 3 — 1. Furthermore, since It is 
finite, there are two unbounded stable regions in rjt . Therefore 

1 2 1 ^ 

E[It+k\It\ < g pi 11 It + 2 p{lt + (j ~ !) l(i + 5^p£(i - M k- i ) 

— — i= 1 i= 1 

< + ^ + If 

( 4 ) 

Our next goal is to show that there exists k for which the constant in front 
of It is smaller than 1. In order to do this, we compute p J k and pf(n,m ) for 
0 < n,m < 2k for some values of k. 

First we take k = 1. Let r] € A. Then P(<r^i( T? )(0) = 0) only depends on r)(x) 
and a v (x) for — 2 < x < 2. So to find the probabilities p ] k and p k (n,m), we 
can just check all possibilities. In the table below we listed them for the case 
a v ( 0) = 0, without loss of generality assuming that 7/(0) = 0 and omitting 
(symmetrically) equivalent cases: 


(»z(*))*=i 2 

K (*))£=- 2 

P(^iw(°) = °) 

00000 

00000 

1/2 

00001 

00001 

1/2 

00010 

00010 

1/2 

00010 

00011 

3/8 

00011 

00011 

5/8 

10001 

10001 

1/2 


It follows that p[ = 5/8 and p[ n = 1/2. As a second example, we compute 
pf( 1,2). So here we maximize P(a R i( r? )(0) = 0) over configurations p for 
which <j J) (—2) = 0 and cr r] (x) = 1 if —1 < x < 2. Again we assume that 
7/(0) = 0. Then the following cases are possible: 


(v{x))%Z 2 -2 

F (^l(d(°) = °) 

01001 

0 

01010 

0 

01011 

0 

10010 

1/2 

10011 

1/2 


Therefore, pf (1,2) = 1/2. For other values of n and m a similar calculation 
leads to the following values of pf(n, in): 



m = 0 

m = 1 

777. = 2 

n = 0 

1/2 

3/4 

1/2 

n = 1 

3/4 

1/2 

1/2 

n = 2 

1/2 

1/2 

0 
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Using ((2D the maximal value of the sum in Q turns out to be equal to 2 
(for g = 4), whence 


E[/ m |/ f ] < ^3 • 2 + 3 ' 8 + 3 ' 2 ) 4 Jt - 


For other values of k, we can proceed analogously. The key observation 
is that P(ir R fc( r; )(0) = 0) is completely determined by rj(x) and a v (x) for 
—2k < x < 2k. For k = 4, we find the desired inequality E[I t+ k\It] < dt 
with c < 1. Below is a summary of the computational results. 


E[/ t+2 |/ t ] < 
nh +3 \it\ < 

Wt+A\It\ < 


(1 29 , 2 61 | 1 191 t _ 121 t 

V 3 ’ 64 3 ’ 128 ' 3 8 ) * 96 

(1 5037 i 2 2687 , 1 2495 1 j _ 55705 j 

\ 3 ' 16384 ' r 3 ' 8192 " r 3 ' 10247 49152-'*’ 

(1 L5371121 , 2 518955 , 1 2371247 1 j _ 200344049 j 

V 3 ' 67108864 ^ 3 ' 2097152 " r 3 ' 1048576/ * 201326592 ±t ' 


The values of pf(n,m) can be written as fractions. Their numerators are 
given in the table below. Entries in the same column have the same denom¬ 
inator, which is written in the second line of the table. All denominators are 
powers of 2. For example pf( 5, 6) = pf(G, 5) = ^|(= 



n =0 

n= 1 

n=2 

n =3 

n=4 

n=5 

n =6 

n=7 

n =8 

denom. 

536870912 

67108864 

67108864 

4194304 

131072 

2048 

512 

16 

i 


(=2 29 ) 

(=2 26 ) 

(=2 26 ) 

(=2 22 ) 

(=2 17 ) 

(=2 n ) 

(=2 9 ) 

(=2 4 ) 

(=2°) 

m =0 

109921252 









m= 1 

125921345 

14557783 








m =2 

116330096 

17164756 

15205719 







m =3 

131398816 

15993296 

18114414 

926907 






m=A 

118712000 

17449960 

14234240 

1036746 

21623 





m =5 

128603456 

15641520 

16272448 

860952 

27214 

307 




m=6 

112440704 

16537568 

12895872 

948720 

17728 

358 

47 



m =7 

119681024 

14745472 

14295552 

773184 

21120 

256 

62 

1 


m =8 

105200384 

14745472 

11496192 

773184 

14336 

256 

32 

1 

0 


Now we have the inequality 


Wt+A\It] < clt, 


with c 


200344049 

201326592 


Taking expectations, 


E [h+A] = EpE[/ t+4 |/ t ]] < cE[I t ], 


Therefore, for all t € N, we obtain 

E[/ 4t ] < c 4 E[/ 0 ] < c\2M + 1). 
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Since ht, is positive integer-valued and the events {ht > 1} are decreasing, 
we get 


P(3 1 : ht = 0) = 


> 


Q{/ 4 , = 0} =l-P f]{/ 4t > 1 } 


\t =o 


\t =o 


1 — lim P ({ht > 1}) > 1 — lim E[/ 4 t] 

t->-oo t—too 

1 - lim c t (2M + 1) = 1, 

£—>• oo 


where we used Markov’s inequality. If ht 0 = 0, then r]t is stable for all 
t > 4to- Therefore 


P(3to : rjt is stable for all t > to) = 1. 


□ 
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